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A quantum error correcting code is a subspace C such that allowed errors acting on any state in C 
can be corrected. A quantum code for which state recovery is only required up to a logical rotation 
within C, can be used for detection of errors, but not for quantum error correction. Such a code 
with stabilizer structure, which we call an “ambiguous stabilizer code” (ASC), can nevertheless be 
useful for the characterization of quantum dynamics (CQD). The use of ASCs can help lower the 
size of CQD probe states used, but at the cost of increased number of operations. 
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I. INTRODUCTION 

In quantum information processing, a quantum error 
correcting code (QECC) is a subspace C, carefully se¬ 
lected to protect from certain noise, any initial state 
1^) € C [l|. Let {Iji)} be a n-qubit basis for C, en¬ 
coding fc-qubit states |j), with 0 < j < 2^ — 1. Such 
a code is a [[u, k]] QECC, where k is the code rate. In 
this work, we will represent the noise using the error basis 
given by elements Ek of the Pauli group P„, the set of all 
possible tensor products of n Pauli operators, with and 
without factors ±1,±L Thus eI = Ek and [EkY = Im 
the identity operator over n qubits. 


A. Stabilizer formalism 

A stabilizer description of error correction is connected 
to classical error correcting codes over GF(4) 0. Where 
applicable, the stabilizer formalism is advantageous in fo¬ 
cussing attention on measurement operators, which can 
be compact, rather than on code words, which can be 
large. A state \ifL) is said to be stabilized by an operator 
S if S'|'0 l) = IV'i)- Let G he a. subset of n — k indepen¬ 
dent, commuting elements from Vn- A [[n, fc]] QECC is 
the 2^-dimensional simultaneous -|-l-eigenspace C of the 
elements of G- A basis for C are the code words IjT)- 
The set of 2”“^ operators generated by G constitute the 
stabilizer S. The centralizer of S is the set of all elements 
of Vn that commute with each member of S: 

z = {PGVn\ysGS,[p,s] = o}, ( 1 ) 

while the normalizer of S is the set of all elements of Pn 
that conjugate the stabilizer to itself: 

Ar= {Pg I =5}. (2) 
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We note that S C Af because the elements of S are uni¬ 
tary and mutually commute. Similarly, Z C J\f because 
elements of the Z are unitary and commute with all ele¬ 
ments of iS. To see that J\f Z, we note that if G Af 
then given any S G S, NSN'^ = S', for some S' G S, or 
NS = S'N. For Pauli operators, NS = ESN, meaning 
S' = ±S. But if S' = -S, then NSN'f = -S, which 
would require that both S and —S are in S. However if 
S G S, then —S is not in the stabilizer, so the only possi¬ 
bility is S' = S, and we obtain that for each N and any S, 
[iV, S'] = 0, i.e., J\f C Z. It thus follows that here Z = Af. 
We have SNljp) = NSIJl) = A^Ijl), which shows that 
the action of N maps code words to code words, and thus 
has the action of a logical Pauli operation on code words. 

A set of operators Ej G Vn constitutes a basis for 
correctable errors if one of the following conditions hold: 

EjEk G S (3^) 

3GgG :[E,Ek,G]^0. (3b) 

The case (I3a1) corresponds to degeneracy. Here 

{'ipL\EjEk\'ipL) = = I, meaning that both errors 

produce the same effect, and the code space is indifferent 
as to which of them happened. Thus either error operator 
can be applied as a recovery operation when one of them 
occurs. The case (l3bl) corresponds to EjEk ^ Af- In that 
case, 3G G G ■ EjEkG = —GEjEk, which ensures that 
G anti-commutes with precisely one of the operators Ej 
and Ek. Thus the noisy logical states Ej\'tpi,) and Ek\'ifL) 
will yield distinct eigenvalues (one being -|-1 and the other 
— 1) when G is measured. The set of n—k eigenvalues ±1 
obtained by measuring the generators G forms the error 
syndrome. The consolidated error correcting condition 
([3]) can be stated as the requirement 

EjEk^Af-S, ( 4 ) 

for every pair of error basis elements, with j ^ k. 

B. Noise characterization and stabilizer codes 

Characterization of the quantum dynamics (CQD) 
of a quantum system is vital for practical applications 
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in quantum information processing, communication and 
computation, on account of environmental decoherence. 
If a represents the quantum state of the system at time 
t = 0, then it evolves under the action of the noise to 

(5) 

m,n 

where x is the process matrix, a Hermitian operator satis¬ 
fying the properties Xi,kE]Ek = I, and Ym Xm,m = 
1 0. The (positive) elements XjJ are probabilities for 
errors Ei to occur. The terms Xj,k (j ^ k) refer to the 
coherence between two distinct errors. 

Quantum process tomography (QPT) denotes a CQD 
technique where for selected input states aj, complete 
state tomographic data Xj^k = Er{aj£(ak)) is obtained. 
The process matrix is derived by inversion of this exper¬ 
imental data. There have been several QPT techniques, 
like standard quantum process tomography (SQPT) [Jja 
and ancilla-assisted process tomography (AAPT) Q. 
Other CQD methods, which bypass state tomography, 
include “direct characterization of quantum dynamics” 
(DCQD) a and “quantum error correction based char¬ 
acterization of dynamics” (QECCD), introduced by the 
present authors (3|. Other developments include an effi¬ 
cient method for estimating diagonal terms of the process 
matrix using twirling Q, which is useful for determining 
QECCs Q. Other related works on channel estimation 
include Ref. 0, a technique like that in Ref. H ex- 
tended to cover off-diagonal Xjk terms, and Refs. |ll|. 

QECCD brings about a twist to the theme of using 
QECCs in that the codes are used not just for protect¬ 
ing quantum states, but also for CQD. This permits one 
to implement CQD concurrently with quantum compu¬ 
tation, making QECCD work “online”. The intuition 
is that the statistics of errors detected during the er¬ 
ror correction process are used to characterize the noise. 
QECCD makes use of the properties of a class of sta¬ 
bilizer codes in which the allowed Pauli error operators 
form a group. 

C. Code ambiguity 

A natural extension of the concept behind QECCD 
would be to adapt stabilizer techniques purely or primar¬ 
ily for CQD, rather than for quantum error correction. 
Freed from correction duty, codes are no longer bound 
to obey (|31). This can be exploited to design codes with 
code lengths smaller than permissible under ([3]), thereby 
making code words easier to implement experimentally. 

However the price to pay for violating m is that some 
errors will be indistinguishable, making stabilizer mea¬ 
surement outcomes ambiguous. The new class of stabi¬ 
lizer codes that we introduce here are therefore called 
ambiguous stabilizer codes (ASCs). In an ASC, the fi¬ 
nal state after recovery may contain a residual logical 
Pauli operation with respect to the initial logical state. 
An ASC generalizes the concept of a degenerate code, 


which is the special case where the only residual logical 
operation after recovery is the trivial one. 

Thus, the purpose of invoking ambiguity-indeed the 
principal motivation behind the construction of ASC’s- 
is to exploit the stabilizer formalism and structure for 
the construction of codes that are better suited for error 
characterization rather than for error correction. 

Code ambiguity entails, as detailed below, that more 
state preparations involving other ASCs are required to 
unambiguously determine the process matrix. Later we 
will find that the number of ASCs required for full charac¬ 
terization scales linearly with ambiguity, in a way made 
precise later. Thus there is a trade-off between spatial 
resources (length of code words) and temporal resources 
(number of operations). We call this ambiguous exten¬ 
sion of QECCD as “quantum ASC-based characteriza¬ 
tion of dynamics” (QASCD). From an experimental view¬ 
point, the above trade-off means QASCD helps simplify 
quantum state preparation at the cost of increased num¬ 
ber of trials and classical post-processing. 

QASCD is, unlike traditional methods of process to¬ 
mography but like the techniques presented in Refs. 
[^ □], a direct method in that it does not require the 
full state tomography of probe states used for CQD. At 
the same time, QASCD may require smaller codes than 
these direct techniques. In particular, to characterize m- 
qubit noise, the above direct techniques require probes 
to be 2m-qubit states or larger. By contrast, with ASCs 
one can beat this bound. For example, to characterize 2- 
qubit noise, one can in principle use (a family of) 3-qubit 
ASCs. 

The remaining article is structured as follows. After 
first developing a theory of ASCs in Section [Hi we study 
in Section imi their specific group theoretic properties as 
would be useful for CQD. In Section IlYl we detail the 
protocol that would be used for CQD by employing (a 
family of) ASCs. The resources, in terms of number of 
ASCs and operations required for CQD, are discussed 
in Section 13 A trade-off between the space resources 
(length of codes used) vs time resources (required con¬ 
figurations) is discussed here. After illustrating our new 
method as applied to a toy 2-qubit noise in Section IVll 
we finally present conclusions in Section IVlII 


II. AMBIGUOUS STABILIZER CODES 
A. Definition and basic features 

A 2^-dimensional subspace C of n qubits, together 
with an allowed set E of Pauli error operators, is ambigu¬ 
ous when two or more errors cannot be distinguished via 
syndrome measurements on the logical state. The indis¬ 
tinguishable errors may require different recovery opera¬ 
tions. Thus ambiguity generalizes the concept of degen¬ 
eracy, and in general prevents error correction. 

An ambiguous set is a collection of allowed Pauli 
errors that cannot be mutually distinguished by a syn- 
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drome measurement. Let = {E^\e'^\ • • • , 

(p) 

where Ey is the jth error in the pth ambiguous set, and 
there are v{p) elements in the pth ambiguous set. Then 
any two errors and where j ^ k, produce the 
same error syndrome. Thus, Pauli errors that are ele¬ 
ments of the same ambiguous set are mutually indistin¬ 
guishable. However, they may not be degenerate, and 
thus may require different recovery operations, making 
an ambiguous code unsuitable for error correction. 

Ambiguity of the code can be represented by par¬ 
titioning E into ambiguous sets. The collection of 
all ambiguous sets is the ambiguous class A = 
with IJpA^P^ = E. The order of 
ambiguity of the code is cr, the number of ambiguous sets 
in the ambiguous class. The degree of ambiguity^ denoted 
7 , is the largest number of Pauli errors in E that can be 
mutually ambiguous, i.e., 

7 = max|A(P)|. ( 6 ) 

p 

A conventional (unambiguous) stabilizer code is charac¬ 
terized by 7 = 1, whereas for an ASC, 7 > 1. Any set of 
up to s known errors drawn from distinct ambiguous sets 
A(p) can be detected, and if the errors are known, they 
can be corrected. 

Within an ambiguous set A^p\ the error elements pro¬ 
duce the same error syndrome. This means that the 

(v) (r>) 

action of two ambiguous errors En' and Em must be 
related by 

Ey^jL)=NEy^jL), (7) 

where A^ is a normalizer element. Note that NE^ \ji,) = 
±E^^^N\jl). Thus, 

2fc_i 2'=-! 

E E^r^^\dL){jL\E^y^ = E Ey)N\jL){jL\NEy) 

j=0 j=0 

2'=-! 

= E 

j'=0 

= n(p\ ( 8 ) 

where j' is simply a re-ordering of j. In other words, 
every element E^ generates the same erroneous sub¬ 
space, with projector However, individual code 

words are not necessarily mapped to the same erroneous 
code word, in view of ©■ Further, from Eq. 0 , we 
have N = E^P'>E^P\ If then 

N' = Elfi>Elf\ Thus, A^t = Af'. 

Note that if [En\E^] = 0, then A^l = A^ (Hermicity 
condition) and thus N = N'. Conversely, if A^ = A^', 
then E^X^E^f^^ = E^r^^E^^\ and thus [E^r!'\E'^^] = 0. 
If {Elf\Em^} = 0, then A^f = —N (anti-Hermiticity) 
and thus N = —N'. Conversely, if A^ = —N'^ then 
E^y^E^y^ = -E^y^E^y\ and thus = 0. 


In contrast to the case with subspaces generated by 
ambiguous errors, projectors to distinct unambiguous er¬ 
roneous subspaces are orthogonal: 

n(p)n(9) = 0, (9) 

li p ^ q. Thus two or more errors belonging to distinct 
ambiguous sets can always be disambiguated. 


B. Ambiguously detectable errors 

Ambiguous errors Em'^ and E^'^ that are linked in Eq. 
© with N = II, where /l is the logical Pauli identity op¬ 
erator, require the same recovery operation. Ambiguous 
errors related by non-trivial logical Pauli operations will 
require distinct recovery operations. Thus, an ambiguous 
code cannot be used for quantum error correction. 

For ASCs, the error correcting condition ([3]) becomes: 

Py^q ^ E^^E^^ i AT, (lOa) 

P = q^ Ey^E'fy e M. (10b) 

Eq. (llOal) implies that quantum error correction can 
be implemented for any collection of known errors which 
belong to distinct ambiguous sets. Eq. (llObI) implies 
that any pair of errors belonging to the same ambigu¬ 
ous set will produce the same syndrome, and thus be 
indistinguishable. In particular, if E^En'^ € S, then 
{i)L\Em''En^ifh) = meaning that the two er¬ 

rors are mutually degenerate, and the ambiguity is harm¬ 
less in the sense that the recovery operation for any one 
of them works for the other, too. On the other hand, if 
Em^En'^ € Af — S, then the erroneous code words they 
produce are related by non-trivial logical Pauli opera¬ 
tions Eq. ©, and the error correcting conditions are 
violated. If one implements a recovery operation favoring 
a single error in each ambiguous set, this will in general 
produce a mixture of states within the code space C', 
which are logical Pauli rotated versions of each other. 

In A, each ambiguous set A^p'> corresponds to the same 
error syndrome, so that order a < 2"“^. By definition, 
the set will contain the element I and, by virtue 
of Eq. (llObI) . only elements of the normalizer Af. The 
remaining sets A^^\A^'^\ ■ ■ ■ will contain Pauli operators 
not present in Af, since they will fail to commute with at 
least one stabilizer generator. 

For unambiguous (and non-degenerate) recovery using 
a linear QECC, the dimension of the code space and the 
volume |E| must satisfy the quantum Hamming bound, 
2'=|E| < 2’^, or 

log(E) < n — k. (11) 

An ambiguous code may violate (fTTl) . though not neces¬ 
sarily. A QECC that saturates Eq. (fTTl) is called perfect. 
The 5-qubit code of Ref. [l^ is such an example. 
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C. Constructing ASC’s 

The simplest way to produce an ASC is by error over¬ 
loading a stabilizer code. This involves allowing addi¬ 
tional errors in violation of condition such that in¬ 
stead condition cni) holds. Ambiguity produced by error 
overloading a perfect code will result in a violation of the 
quantum Hamming bound (11111 . while for an imperfect 
QECC, a sufficiently large amount of error overloading 
would be required to violate Ineq. dni). For example, 
consider the (perfect) 5-qubit code of Ref. 

|0l)5 = (-100000) -H 10111) - 110011) -f IlllOO) 

-HlOOllO) -H 101001) -f IIOIOI) -h 111010)) 

|li)5 = j— (-111111) -H IIOOOO) -h 101100) - 100011) 
2 v 2 

-h|11001) -t- IIOIIO) - 101010) - 100101)), (12) 

which corrects an arbitrary single-qubit error on any 
qubit. The code space is stabilized by generators 
IXXYY,IYYXX,XIYZY,YXYIZ. They can each 
take values ±1, thereby determining 16 syndromes, cor¬ 
responding to the 16 allowed errors E = {I,Xi,Yi, Zi} 
where i = 1, • • • , 5. By allowing any more errors into the 
error set E, we introduce ambiguity, and also violate (Ell¬ 
in Table HI we present a partial listing of the ambiguous 
class A for this code. In all, it has 1 + (5) -3-1- ( 2 ) -3^ = 106 
arbitrary 1-qubit and 2-qubit errors, of which 49 are dis¬ 
played. The errors are partitioned into their respective 
ambiguous sets, labelled by the corresponding error syn¬ 
drome. Set has only 1 element, I, since all other 
elements of Af have a Hamming weight greater than 2. 


Another way to create an ASC from a stabilizer code 
is by syndrome coarse-graining: dropping one or more 
syndrome measurements. For example consider not to 
measure the last stabilizer of the code ED- From the 
first entry corresponding to syndromes (the wn-error- 
overloaded case) of the TableUit can be seen that |.A| = 8 , 
= {/, Ai} corresponding to syndrome (-1- -I- -I-), 
= {Yi,Zi} corresponding to syndrome (-I--I —), and 
so on. The order of ambiguity is halved and the degree 
of ambiguity is doubled. 

A final method to obtain an ASC begins by construct¬ 
ing a stabilizer code that corrects arbitrary errors on m 
known coordinates. An ASC may then be obtained by 
allowing noise to act on m' known coordinates, where 
m' > m. A detailed description of this method and its 
application to the characterization of quantum dynamics 
H are considered below. 


III. AMBIGUOUS GROUP 

An arbitrary error on I qubits can be expressed as a 
linear combination of 4* Pauli operators. Suppose these l- 


-h-h-h-h 

-h-h-h - 

-h-h- -h 

-h-h- 

I 

Xl 

Yl 

Zl 


Y2Y3 

Z2Z3 

X2X3 


X3Y4 

Y3X4 

Z3Z4 

-h—h-h 

+-- h- 

H-h 

+ - 

A2 

Ys 

Y4 

X3 

^1X3 

X1X2 

Y1X2 

^1X2 

Y3Z4 

Z2Y3 

Y2Z3 

Z2X4 

— h-h-h 

-++ - 

-- h- -h 

-- h — 

Y3 

Y2 

X4 

X5 

A1Y2 

X1Y3 

Z1Y2 

Y1Y2 

A 2^4 

^ 3 Y 4 

Z2X3 

X2Z3 

- h-f 

--h- 

- h 

— 

^4 

Z2 

^5 

^3 

A1Z2 

Y1Z3 

Z1Z2 

Yi^2 

A2Y3 

X3X4 

X1Z3 

Y2Y4 


TABLE I. Ambiguous class (partial listing) for the ASC ob¬ 
tained by error-overloading the code E, to allow arbitrary 
errors on any two qubits. Each error syndrome labels an 
ambiguous set. The first error row in each column corre¬ 
sponds to arbitrary single-qubit errors allowed in the original 
QECC. Inclusion of the two-qubit errors (second and third 
rows of the table) to the list turns the QECC into an ASC. 
In all, there are 106 elements in the ambiguous class, with 
I I = 1 and | | = 7 for p = 1, 2, ■ • ■ , 15. Thus the degree 

of ambiguity is 7. For example, the full ambiguous set, cor¬ 
responding to the syndrome -I- -|- H— has four more elements 
= X4X5,Ei^'> = Z 3 Z 5 ,eP =X2Y5, and E^^^ = Z 2 Z 4 . 
The normalizers between Eq^^ = Xi and other elements in 
the set are XYYII Zl, XIXYI -Yl, XIIXX Zl, 
XIZIZ -Xl, XXIIY -Yl and XZIZI -Xl. 

Any set of sixteen elements, with one drawn from each am¬ 
biguous set will satisfy condition Eq. (llOall . while any pair 
of errors within a column satisfy Eq. (llObI) and thus are am¬ 
biguous. Further note that the product of ambiguous errors 
linked by the same logical Pauli are mutually degenerate (e.g., 
E^^'^ E^'’ € S), and are correctable by the same recovery op¬ 
eration, while those linked by different logical Pauli operators 
are not (e.g., E^''^E^^'' £ M — S). 


qubits form a subsystem of n qubits prepared in a [[n, /c]] 
stabilizer code. Setting |E| := 4^ in Ineq. (fTTl) we find: 



(13) 


This means that a 5-qubit code can correct all possible 
errors on at most 2 fixed coordinates. An example of a 
perfect code of this kind will be presented later. We thus 
obtain a [[n, fc]] ASC by allowing m noisy coordinates, 
where m > / in Ineq. m- The order a of the code is 
just the number of syndromes, 2 "“^, while the degree of 
ambiguity 7 = 4™/2’"-'= = 22 '"""+'=. 

Suppose we are given a [[n, fc]] ASC with errors allowed 
on m known coordinates. It is worth noting here that the 
set of errors (including the factors ± 1 , ±i) forms a group, 
i.e., E = Vru- Furthermore, the subset 25 of Vm that is 
ambiguous with (the trivial error on the m qubits) 
constitutes a group, the ambiguous group, as shown be¬ 
low. 
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Theorem 1 Given a [[n, fc]] ASC with E = Vm, the sub¬ 
set 25 of allowed errors that correspond to the no-error 
syndrome forms a normal subgroup. 

Proof. Note that if Bj,Bk € then Bj and Bk both 
commute with all stabilizers, by virtue of Eq. (jlObI) . 
(Note that this doesn’t imply that [Bj, Bk] = 0. Thus the 
subgroup is not Abelian.) For any element G G G, then 
[BjBk,G] = BjBkG—GBjBk = 0, meaning that BjBk € 
25. This guarantees closure of the set. By definition, 
is an element of this set, and a Pauli operator is its own 
inverse. Thus all required group properties are satisfied. 
Normalcy of the subgroup (the equality of the left and 
right cosets) is guaranteed because we implicitly include 
Pauli operators with and without factors ±1 and ±b ■ 

For an ASC obtained in this way, the ambiguous class 
A has a simple structure. It corresponds to a partition 
of Vm, determined by the quotient or factor group 

(w) 

This means that any element B in Vm is either in 25 or 
can be expressed as the product of an element in *8 and 
an element not in *8. 


A. Example of a [[3, 1]] ASC 

A [[3,1]] perfect QECC that unambiguously corrects 
errors on the first qubit is Q: 

|0L)3 = i(|001) + |010) + |100) + |lll)) 

|1l) 3 = i(|110) - |101) + 1011) - 1000)), (15) 

whose stabilizer generators are given by the set Ga = 
{XIX,YYZ}. The stabilizer is thus the set of four ele¬ 
ments, S 3 = A X 2^ = A X {/, XIX, YYZ, ZYY}, where 
the pre-factor indicates possible factors ±1, . The nor- 

malizer A 3 is the set of all elements of V 3 that commute 
with the elements of 53 . (We note that a Pauli operator 
P commutes with every element of 1 S 3 iff P commutes 
with every element of Ga-) 

For code (HH), the normalizer A 3 is given in Table |lll 
Since there are only four logical Pauli operator, various 
normalizer elements map to the same logical Pauli oper¬ 
ator by virtue of their effect on the code words (fT51) . The 
subset 53 (the first column) corresponds to the identity 
logical Pauli operation II, while the elements of A 3 — S 3 
correspond to non-trivial logical Pauli operations, as tab¬ 
ulated in the remaining columns of Table HU 

We create an ASC for the code (fTSll by allowing errors, 
in addition to the first coordinate, also on the second 
coordinate. There are four elements in Table HD that have 
no non-trivial operator on the last qubit, i.e., they are 
elements of P 2 ® II 3 , where I3 is the identity operator on 
the third qubit. They are {III,XZI,IYI,XXI}, which 


Il 

-Xl 

Yl 

Zl 

III 

XZI 

IYI 

XXI 

XIX 

IZX 

XYX 

IXX 

YYZ 

YXY 

YIZ 

YZY 

ZYY 

zxz 

ZIY 

ZZZ 


TABLE II. Normalizer for the [[3,1]] stabilizer code (1151) and 
the logical Pauli operations they map to. All elements com¬ 
mute with the elements of 53, while their commutation prop¬ 
erties amongst themselves reflect the logical operation they 
map to. Thus, an element in the column Yl commutes with 
all elements in the same column and those in the column II, 
but will anti-commute with every element in the columns -X^ 
and Zl- On the other hand, the elements in the column II, 
which are precisely those of 53, commute with every other 
element in the normalizer. 


-b-b 

-b - 

- -b 

— 

Normalizer 

I 

Xl 

Yl 

^1 

II 

Y2 

X1Y2 

Y1Y2 

ZiTz 

Yl 

X1X2 

X2 

Z1X2 

Y1X2 

Zl 

X1Z2 

Z2 

Z1Z2 

Y1Z2 

-Xl 


TABLE III. Ambiguous class A 3 for errors on the first 2 qubits 
of 3-qubit code (1151) . depicting the quotient group (1161) . The 
first column is the ambiguous group QSs, drawn from Table 
im subject to the requirement that the operator on the third 
qubit is identity. The remaining three columns are its cosets 
AilBs, Yi ^3 and ZitBs, which represent ambiguous sets. The 
last column lists the normalizer element with respect to first 
element in the column, in the sense of Eq. 0. For exam¬ 
ple, the normalizer element that maps error Z1Z2 to error 
Yl or vice versa is Nz^z^.Yi = Nyi,ZiZ 2 = —Xl, while the 
normalizer element which maps error yiy 2 to error Z 1 X 2 is 
XziX 2 ,YiY 2 = ZlYl = iXl, while Nyj^Y 2 ,ZiX 2 = YlZl = 
-iXk. 


constitute the ambiguous group 253 . The partitioning of 
A for ASC (fT^ with £ = 7^2 can be represented by the 
quotient group: 


This is depicted in Table uni where the first column is 
the ambiguous subgroup * 83 , and the other columns are 
its cosets and the other ambiguous errors. 


IV. APPLICATION TO NOISE 
CHARACTERIZATION 

While code ambiguity makes ASCs not useful for quan¬ 
tum error correction, they can be used for experimentally 
studying noise. In this Section, we elaborate on the intu¬ 
ition presented earlier, of extending QECCD by replacing 
the use of stabilizer codes with that of ASCs. 
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A. Ambiguity and QEC channel-state isomorphism 


The basis of QECCD is the quantum error correc¬ 
tion (QEC) isomorphism, qualitatively similar to the 
Choi-Jamiokowski channel-state isomorphism, which as¬ 
sociates a correctable noise channel with the unique erro¬ 
neous logical state corresponding to a given input logical 
state. This clearly is necessary if complete information 
about the channel is to be extracted via measurements. 
In the presence of ambiguity, for any initial logical state, 
it can be shown that one can always construct two or 
more noise channels such that they produce the same er¬ 
roneous logical state. Thus QEC isomorphism no longer 
holds. 

In QECCD, the QEC isomorphism is leveraged 
through some state manipulations to yield full noise data. 
The basic idea is that the syndrome obtained from the 
stabilizer measurement is used to correct the noisy state, 
while the experimental probabilities of syndromes will 
characterize the noisy quantum channel. While direct 
syndrome measurements yield the diagonal terms of the 
process matrix, for off-diagonal terms preprocessing via 
suitable unitaries is required. For the purpose of noise 
characterization, the code qubits are divided into two 
parts; (a) the qubits on which the elements of E act non- 
trivially; (b) the remaining qubits. 

The former qubits constitute the principal system P, 
whose unknown dynamics is to be determined. The latter 
qubits constitute the CQD ancilla A, and are assumed 
to be clean, i.e., noiseless. Suppose the full system P-t-A 
is in the state 


2'“-! 

j=o 

where {Ijl)} denotes a logical basis for the code space 
of a [[p -I- q, fc]] ASC (which encodes k qubits into n = 
p + q qubits) such that allowed errors in the p known 
coordinates of P can be ambiguously detected. 

The main difference QASCD has with respect to 
QECCD is that QASCD employs more than one code to 
fully characterize the noise. Herebelow, we present de¬ 
tails of the QASCD protocol, which has a quantum part, 
which is experimental, and a classical part, which con¬ 
cerns post-processing data from experiments. The quan¬ 
tum part involves using state preparations and syndrome 
measurements of different ASCs to determine Xm,n am¬ 
biguously. The classical part involves simultaneous equa¬ 
tions to disambiguate the ambiguous experimental prob¬ 
abilities. 


B. Direct measurement 

Let Q be an ASC that can detect noise £, with asso¬ 
ciated process matrix y. Let Ea^ (j = 0,1, 2, • • • , 7 — I) 


be the elements of an ambiguous set in Q, with de¬ 
noting any one of these a^’s. It is convenient to employ 


the notation 



EoUl)- The probability that one 


of these ambiguous errors occur: 


M = Tr f (|'0 l)('0l|) 


2^-1 




3=0 


2 "-! 




3=0 

+ + • • •] \JI) 

= • • • + Xa,.c, + Xa,.a. (V'A 


\0i2^0i2 


= -HXai.ci +Xai.a2('0L|^l,2|V’L) + Xas.ci (V'L 1^2,1 IV'i) 

H” \ciL2-,OL2 “I” ' ’ ’ 

~ 'y y Xaj.aj + 2 Re (xaj.afc (Ajjfe)/,) , (18) 




where Nm,n = E^En- Note that because Em and 
produce the same syndrome by virtue of ambiguity, Nm,n 
so defined will commute with all elements of the stabi¬ 
lizer. 


Let D = 2P, the dimension of P. In an unambiguous 
code, the diagonal terms of y would appear as prob¬ 
abilities of syndrome measurements [^. Now, however, 
any measurement outcome probability will contain con¬ 
tributions from the probabilities of 7 ambiguous errors 
plus ( 2 ) off-diagonal terms between these ambiguous er¬ 
rors. Of the 4^ can be disambiguated by using as many 
different initial state preparations, by exploiting the fact 
that the y terms have factors given by expectation val¬ 
ues of different normalize!' elements (logical Pauli oper¬ 
ations). However, the problem of disambiguation would 
still remain within each such ‘logical Pauli class’, i.e., dif¬ 
ferent pairs of ambiguous errors {Ej,Ek) such that the 
normalizers EjE^ correspond to the same logical Pauli 
operation. This is related to limits imposed by the am¬ 
biguity and can be sorted out by using other ASCs. For 
accessing cross terms terms of y across ambiguous sets, 
we use the unitary pre-processing described below in Sec¬ 
tion IIV Cl and IIVDI based on the method introduced by 
us in Ref. Q. 

As an example of result (fT51) . for the data in Table HTTl 
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the probabilities to obtain all the outcomes are 

p(++) = Xij + XF 2 .Y 2 + XxiX2,XiX2 + XxiZ2,XiZ2 

+ 2 X [(Re(x/.XiZ2) + ^™-iXY2,XiX2)){XL) 

- (R-e(x/.y-2) - ^^{xxiX2,XiZ2)){Yl) 

+ (Ini(x>2,XiZ2) + ^(^{xi,XiX2)){Zl)], 

-) = XXi,Xi + XXiY 2 ,XiY 2 + XX2,X2 + XZ 2 ,Z 2 

+ 2 X [(Re(xxi,Z2) - Im(xxiy2,^2))(^i) 

+ (R-e(xxi.XiX2) - Ini(xx2.Z2))(^i) 

+ (Ini(xxiyi,Z2) + R'e(xxi,X2))(-^i)], 

pi —1“) = XXi.Xi + XXiX2,XiX2 + XZiX2.ZiX2 + XZiZ2,ZiZ2 

+ 2 X [(Re(xxi.ZiZ2) + Iin(xxiX2,ZiX2))(-’^L) 

+ (Ini(xziX2,ZiZ2) - Re(xxi,XiX2))(^i) 

+ (Im(xFiy2,ZiZ2) - R-e(xxi.ZiX2))(^L)], 

pi -) = XZi,Zi + XZiY2,ZiY2 + XFiX2,XiX2 + XYiZ2,YiZ2 

+ 2 X [(Im(xziF2.XiX2) - Re(xzi.yi22)(-’^i) 

+ (Re(xzi,ZiY2) - Iin(XYiX2.YiZ2)(^i) 

+ +(Im(xziY2,YiZ2)Re(xzi,YiX2))(^L)]- (19) 

By choosing input |0)i, one finds p(++) = xij + 
XY2,Y2 + XX 1 X 2 .X 1 X 2 + XX 1 Z 2 .X 1 Z 2 + 2Re(x7.XiX2) + 

2Im(xY2.XiZ2) = C + 2Re(x7,XiX2) + 2Im(xY2,XiZ2)- 
By choosing input |+)l = + |1)l), one finds 

Pi++) = C + 2Re(x7,XiZ2) + 2Im(xY2.XiX2)- By choos¬ 
ing input |t)L = -^i\0)L -H7|1)l), one finds Pi++) = 

C + 2Re(x7,Y2) - 2Im(xxiX2,XiZ2)- We thus have four 
unknowns, given by C (the diagonal contributions), and 
the coefficients of (Xl), {Yl) and (Zl)- One more in¬ 
put, say cos(0)|Oi) -I- sin(0)|li) will suffice to determine 
these 4 quantities. It will thus suffice to determine C. 
More generally, 4^ (the number of logical Pauli opera¬ 
tions) preparations are needed to solve for C. When 
C is extracted for each outcome, then each code gives 
D'^lx = 2"“^ equations. Note that we have ignored the 
off-diagonal terms for ambiguous errors, since they will 
be dealt with in other ASCs, where they correspond to 
off-diagonal terms that are unambiguous. 

C. Preprocessing with U 

For a given ASC, to derive off-diagonal terms between 
errors in different ambiguous sets, we preprocess the sys¬ 
tem by applying a suitable unitary [/, based on the idea 
we introduced in Ref. Q. However, even this may allow 
one to access only the real or imaginary part of these 
terms. To access the other part, one uses a further pre¬ 
processing described in the next Subsection. 

The unitary U will be in one of two forms. In the first 
form, U = ■^{Ea + EiD, in case [Ea,Efl\ ^ 0. In the 
second form, [/ = -^ (Aa -|-zAb), in case = 0. 

We require Ea and i?t, to be mutually unambiguous in 
the given ASC because otherwise, as explained later, we 
obtain a situation similar to not using U, as far as noise 
characterization is concerned. 


Let us consider the first case. Suppose the pre- 
processed noisy logical state produces an ambiguous out¬ 
come Ej. Let gAjEj = EaE^j, where the E^.'s consti¬ 
tute an ambiguous set, and gA^ G {ilji*} is the Pauli 
factor. Similarly, let gB^Ej = EhEp^, where the EpPs 
constitute an ambiguous set, and gsj G {±1,±7} is a 
Pauli factor. 

When C/(a, b) is applied to the noisy logical state, and 
an outcome x has been observed, then one of the Ej 
must have been detected, and thus the only contributing 
terms of £1 (pT,) will be those restricted to | iplf) and | ). 

Denoting by He the projector to the code space C of the 
ASC, the probability to observe x when U[a, h) has been 
applied is: 

e(a, h, x) = Tr ([/ [£{\^l){^l\)] iE^^cE,)) ■ (20) 

The terms within the square bracket in Eq. dsni) that 
would make a contribution to the probability of obtaining 
ambiguous outcome E^ are: 

• • • + Xai,ai|V'i“'^)(V'i“'^l + Xai,a2lV'i“'^)(V'i“'^l 
+ Xa2.«ilV'i“'^)(V'i“'^l + Xa2.a2lV'i“'^)(V'i“'^l +••• 

+ X«„/3.+ Xft.aJV'f +••• 

+ ••• ( 21 ) 

When the expression in Eq. m is left- and right- 
multiplied by C/(a, 6), then the only resulting terms that 
contribute to the Ihs of Eq. dini) are: 

••• + Xc<i,ai|V'L^)(V’i^l + Xai,a25Ai3A2l^L^)(V'L^I 

+ Xc. 2 ,c^gA 2 gAML^){i’L'’\ +Xc. 2 ,a 2 \^L^){^L^\ + ••• 

+ Xai,/3i5Xi5BjlAL^)(^i^| +X/3i.ai573iflAjV’L^)(^L^I + ••• 

+ Xai,p29A2g*Bj'tpL^){^pL^ \ +X/32.ai57323xJV'L^)(V’L^I 

+ ••• ( 22 ) 

The contribution of the first term in Eq. (na) to the 
probability in Eq. (1^ would be: 

eci.ci = Xai.ai 

i=i 

= Xai ,ai I (23) 

since the traced out quantity has support only in the 
erroneous code space E^C' (i.e., the code space C' shifted 
by the ambiguous error). Analogously, the contribution 
of the fourth term in Eq. (1^ to Eq. (EiH) would be 
ect 2 ,Q 2 = Xa 2 ,a 2 - III fashion, the contribution of the 
fifth and sixth terms in Eq. (0^ to Eq. (OUl) would be 
£qi,/3i = Xcti,/3i Sud epi,c(i = XPl.ai- 

The contribution of the second term in Eq. to the 


probability in Eq. (l20t would be: 

2 '= 

eai.os = Xai,a29A^gX Ol''Vi ^^ I 

i=i 

= Xai.a2 5Ai5A2(V'L 

= Xai.a25Ai5i2(^2.l)L, (24) 

where N 2 ^i is the normalizer element that propagates 
error Ea-^^ on a logical ket to Ea 2 - The contribution 
of the third term in Eq. to the probability in 

Eq. (EiH) would be, analogously to Eq. (IMl) . namely, 
£ 02 ,ai = Xa 2 ,aigA 2 9 2 )L- In like fashion, the con¬ 
tribution of the seventh and eighth terms in Eq. (l22l) 
to Eq. (gni) would be eai,/32 = Xat,p29Ai9B2^Ni^2 )l and 

£/32.ai = X/32,ai9B29Ai{^2,l)L- 

Putting together all these etc., terms into 

Eq. dsni), we obtain: 

^{a,b,x) = i {^Xc.,,c.,+XPBfi2 j 

T ^ ^ [He (Xqj(-^ j.fc)i) 
j<k 

+ ^^{X/3j,Plc9Bj9BkiNj,k)L)] 

+ ^ [Re {x<.,A9A,9ki^Bk)L)] , (25) 

where the (N) terms, being always real, can be removed 
out of the argument of Re or Im. 

In constructing U{a,b), the errors Ea and Et should 
not be mutually ambiguous. Otherwise, the result is ef¬ 
fectively the same as that direct measurement without 
preprocessing using U{a,b). To see this, consider an ap¬ 
plication of this method to Eq. (1191) , with U{X 1 X 2 ,¥ 2 ) = 
-j^{XiX 2 + Y 2 ). We find: 



aj,gA 

Pk,gB 

I 

X 1 X 2 ,1 

Y2a 

¥2 

X 1 Z 2 , i 

7,1 

XIX 2 

7,1 

XiZ2,i 

XIZ 2 

¥2,i 

XlX2,-l 


Erom (l25]l . it follows that with pre-processing by 
U{XiX 2 , ¥ 2 ), the probability expressions in the example 
(O are altered altered, e.g., 

p(++) = 2 X [Re(xxiJf2.V2) + Re(xi,XiZ2)](^L} 

+ 2 X [Re(x/,V2) -I- -Im(xxiX2.XiZ2)](^L) 

-I- 2 X [Re(x/,XiX2) + Ini(xv2.XiZ2)](TL) 

- 2 X [Im(xF2.XiX2) + Im(x/.XiZ2)]- (27) 

Thus, the ambiguous errors in the coefficients of normal¬ 
izer expectation values remain the same even though the 
particular normalizer element changes. 

Now, let U = -^(^1 + ^i): where Xi and Zi are seen 
to be unambiguous for code (ITS!) . Set the outcome to be 


‘-I--I-’. This fixes Ej. Thus: 




•) 9 A 


Ep,gB 


N 

1 

1 

5 

Xi,l 

Xi¥2,1 

772,1 

^ 2,1 

1 

ZiA 

^ 1 ^ 2,1 

yiX2,-i 

¥iZ2,-i 

1 

II 

Yl 

Zl 

-Xl 


(28) 

The coefficient to ^(Xi,Zi,++) can be read off 

(PS)). using (E3, by forming cross-terms between ele¬ 
ments of the second and third columns, such that their 
corresponding logical Pauli operators multiply to Xl up 
to a sign gA- In the present case, this is seen to be 

(Xb) [Im(xxi.ViZ2 - XXir2,r:_X2) + Rkxx2,ZiY2 + XZ2,Zi 

(29) 

We can thus form cross-terms between all ambiguous sets 
using suitable U. 

In the second case, [Ea^Eb] = 0 and we set U = 

^ ^ result, instead of Eq. (1^ . one gets: 

••• + Xai,ai|'0L^)(V'L^I + Xai,a25AiSA2l^L^)(^L^I 
+ Xa2,aiS'A25(4il^L^)(V'L^I +Xa2.c<2l^L^)(V'L^I + --- 
-*Xai./3i5Ai5Bj^/'i^)(^L^| +*X/3i.aigBi5AjV'L^)(V'i^l H 

- *Xai ,/325Ai5 bJ^^ I + *X/32.ai5S25 Ai I ^) (V-i ^ I 

H- 

Consequently, one obtains in place of Eq. (123): 

^(a, b,x) = ^ Xa,,a, + Xft.ft + Xa.-.ft + Xft.a^ ^ 

+ {XaB»k9A^g*AkNjx)L) 

j<k 

+ Re(x/32./3fc5Bi5k(%fc)i)] 

+ {Xaj,gk9Aj9k(^jkL)] , (31) 

t/fc 

where, like before, the (N) terms, which are always real, 
can be removed out of the argument of Re or Im. It is 
worth noting that in Eqs. (123) or dSD, in the terms 
that contain Pauli factors, the matter of whether the 
real or imaginary part of the process element of the pro¬ 
cess matrix contributes to the measured probability, de¬ 
pends on whether the Pauli factors are of same type 
(real / imaginary). 


D. Toggling 

The method of Section llV Cl gives only the real or imag- 
ina^ parts of the cross-terms. Using an idea we proposed 
in [3|, we solve this problem by pre-processing the noisy 
state even before applying U. Consider a density oper¬ 
ator I^a) subjected to the phase operator 
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given by the diagonal T = e*®“|0)(0| + e*®i|l)(l|. Then, 

if 6 q = —9i = 7 r/ 4 , one finds TpT^ = 

meaning that the imaginary and real parts of the off- 
diagonal terms have been interchanged or ‘toggled’ (apart 
from a possible sign change). 

Similarly, now we construct 

(T—1 

(32) 

m—0 

where a is the number of ambiguous sets (order of am¬ 
biguity), n™ is the projector to the erroneous logical 
space given by EmC {Em being any one error from each 
ambiguous set), and 9m S {±f}, with equal entries 
with both signs. Prior to U, we apply the operation 
r+ = T 0 I', where I' acts trivially outside the cor¬ 
rectable space, i.e., the code space C plus the erroneous 
code spaces. 

For example, suppose we construct the toggler T'^ hav¬ 
ing 9e^ = —9ep = ±f, then in place of (^5)) we have: 

{Xl) [Re{xxuYiZ2 - XXiY2,YiX2) + Im(xx2.Zi>2 + XZ2,Zi)] , 

(33) 

i.e., cross-term where p and v come, respectively, 

from ambiguous set E^ and Ep, get their real and imag¬ 
inary parts toggled. 

The tools described in this and the preceding two sub¬ 
sections, as well as the different ASCs, form our reper¬ 
toire for characterizing the noise in the method of ASCs. 


fa ib \ 
-ib* 1 - a J ’ 


V. RESOURCES 

We may begin by supposing that data from 7 ASCs 
will suffice, giving the required equations to solve for 
the variables. These equations will correspond 
to an adjacency matrix, wherein the U ^/7 rows corre¬ 
sponding to each code will sum to a unit row, i.e., one 
with I’s in all columns. Thus there are (at least) 7 — 1 
constraints among the equations. Adding one more 
code will introduce jx equations and one more con¬ 
straint i.e., 2"“^ — 1 constraints. If there are no other 
constraints in the first rows, and if 2 ”“^ — 1 > 7 — 1 , 
i.e., n — k > p, then the remaining required linearly in¬ 
dependent equations can be found from the last code. 
Thus, in general, with 7 -fold full degeneracy, the neces¬ 
sary number of preparations is 7 + 1 . 

More generally, because of the failure of QEC iso¬ 
morphism with ambiguous codes, of the O (4"* x 4"^) 

terms in the process matrix, only O x inde¬ 
pendent terms can be determined per ASC, implying 
that a full characterization would require p = 0 ( 7 ^) dif¬ 
ferent ASC’s. Also, syndrome measurements on each 
ASC yields D'^jx = outcomes. We may thus 

estimate that the number of conhgurations required is 

c = O = O per ASC. Thus in all, count¬ 

ing each ASC as a separate configuration, we require /i x c 


conhgurations, i.e., 0 ( 74 "*), meaning that there is a fac¬ 
tor 7 excess when using ambiguous codes. (Moreover 
each code would require up to 4^ state preparations for 
disambiguation of the Pauli logical classes.) This can be 
considered as a time cost to pay for the saving in ‘space’, 
i.e., in terms of number of entangled qubits used. 

Now we present an example of applying QASCD, with 
three 4-qubit ASCs being used to characterize a 2-qubit 
noise. 


VI. ILLUSTRATION USING A FAMILY OF 
THREE 4-QUBIT AMBIGUOUS CODES 

Consider the [[4,1]] ASC for arbitrary errors on the 
hrst two qubits, constructed by dropping the last qubit 
of the [[5,1]] QECC of Ref. [T||: 

|0i)4 = ^ (-10000) + 10010) + 10101) + 10111) 

-| 1001 ) + | 1011 ) + | 1100 ) + | 1110 )) 

|li)4 = ^ (-11111) + jllOl) + 11010) + 11000) 

-|0110) + 10100) + 10011) + 10001)), (34) 


whose stabilizer generators are XIIX,YIXY and 
YYZZ. The following equation presents two other such 
codes and which are two fold amiguous: 

|0i)4 = Tr|^|0i), |ii)4 = Fi|^|u)4, 

|0i)4 = H^i\QL)A, |li)4 = H^\\Il)a, (35) 

where iJzv = ^(|0)(0|+*|0)(1|+*|1)(0| + |1)(1|), ilyx = 
i((l+*)|0)(0| + (l + *)|0)(l|-(l-*)|l)(0| + (l-*)|l)(l|). 
The corresponding stabilizer generators and the error 
syndromes are given in Table HVl 

By method described in Sec. IIVI the statistics of syn¬ 
drome outcomes on QECs C^, and C^, can com¬ 
pletely determine the process matrix Xm,n corresponding 
to an arbitrary 2-qubit noise £. It can be noticed from 
the Table lIVi that the normalizer corresponds to logi¬ 
cal Yl. By the direct measurement, as in Eq. (HI , we 
get Xa^,a^ +Xa 2 .a 2 + 2Re (x„i,^ 2 ((li))■ By choosing any 
state of other than | f); (Xl) vanishes. The direct 
syndrome measurements on suitably prepared yields 
the following expressions, 

XI,I + XTiYs = “1^ XXi.Xi + XX^Y2,X^Y2 = bi, 

XX2,X2 +XZ2.Z2 = ci,xxiX2,XiX2 + XXiZ2,XiZ2 = di, 
XYiX2,YiX2 + XYiZ2,YiZ2 = ei, XYi,Yi + XYiY2,YiY2 = /l, 
XZiZ2,ZiZ2 + XZiZ2,ZiX2 = 9i,XZi,Zi + XZiY2,ZiY2 = ^i- 

(36) 

Similarly procedure followed on yields 

Xi,i + XZ2,Z2 = a2, XXi,Xi + XXiZ2,XiZ2 = ^ 2 , 

XYi,Yi + XYiZi,YiZ2 = C2, XZi,Zi + XZiZ2,ZiZ2 = d2- 

(37) 




10 



II 

Xi 

X2 

Yi 

Zi 

XX 

YX 

ZX 


Y2 

XY 

Z2 

YY 

ZY 

xz 

YZ 

ZZ 

XIIX 

+ 

+ 

+ 

- 

- 

+ 

- 

- 

YIXY 

+ 

- 

+ 

+ 

- 

- 

+ 

- 

YYZZ 

+ 

- 

- 

+ 

- 

+ 

- 

+ 


II 

Xi 

X2 

Yi 

Zi 

XX 

YX 

zx 


Z2 

xz 

Y2 

YZ 

zz 

XY 

YY 

ZY 

IZZX 

+ 

+ 

- 

+ 

+ 

- 

- 

- 

XIIX 

+ 

+ 

+ 

- 

- 

+ 

- 

- 

YZYZ 

+ 

- 

- 

+ 

- 

+ 

- 

+ 

c 

II 

Xi 

Yi 

Y2 

Zi 

XY 

YY 

ZY 


X2 

XX 

YX 

Z2 

zx 

XZ 

YZ 

ZZ 

IXXZ 

+ 

+ 

+ 

- 

+ 

- 

- 

- 

XIXZ 

+ 

+ 

- 

+ 

- 

+ 

- 

- 

YXYX 

+ 

- 

+ 

- 

- 

+ 

- 

+ 


TABLE IV. Ambiguous class for the three 4-qubit codes. 
The Hadamard operation Hzy [Hyx) toggles errors Z and 
Y (errors Y and X) while keeping error X {Z) fixed, and the 
above syndromes are corresponding toggled versions of each 
other. 

From we obtain the following expressions 

XIJ + XX2,X2 = 03 , XXi,Xi + XXiX2,XiX2 = h, 

XYi,Yi + XYiX2,YiX2 = C 3 , XZi,Zi + XZiX2,ZiX2 = d-3. 

(38) 

The above 16 expressions suffice to determine the di¬ 
agonal terms of the process matrix. To demonstrate how 
the method works for off-diagonal terms, we consider its 
application to the noise 

^a{,Pl)= SpL H-=— {XiPlXi + XZpi^XZ Y 2 PLY 2 

5 

+X2PLX2 + XXplXX) + i ((a + ib)XlpLX 2 

6 

+{c + id)pLXX + {e + if)XZpLY 2 + c.c) (39) 

In the present case, for solving the off-diagonal terms 
using Eq. (ED), the following set of linearly independent 
equations for off-diagonal terms are obtained by perform¬ 
ing unitary operations U(a, b) followed by syndrome mea¬ 
surements on C^, and respectively 

i{I,XiX2,I) = xi,i + XXiX2,XiX2 + XF2.F2 + XX1Z2.X1Z2 

+lm{I,XiX 2 )+Re{Y 2 ,XiZ 2 ), 

^{I,XiX 2 ,Xi) =xxi,Xi +XX 2 .X 2 +XF 2 .F 2 +XX 1 Z 2 .X 1 Z 2 

+Im(Xi,X2)-Re(r2,Xi^2), 

^(/, X 1 X 2 ,1) = XIJ + XJt lX 2 ,XlX 2 + XXi.Xi + XX 2 ,X 2 
+lm{I,XiX 2 )-lm{Xi,X 2 ). 

(40) 

In Eq. (l40l) . the diagonal terms are obtained without 
pre-processing with unitaries using in Eq. ED, dSD and 
(1551) . Solving the above set of equations we obtain the 
off-diagonal terms of the process matrix corresponding to 


Im(/, X1X2) = 1(01+02 + 03) = 

2 o 

Im(Xi,X2) = ^(Oi + O2 - O3) = 

Re(Y2, X1Z2) =0i-^i0i+02 + O3) = ( 41 ) 

where Oi = ^(I,XiX2,I)-{xi,i+XXiX2,XiX2+XY2,Y2 + 
XXiZ 2 ,XiZ 2 )^ C>2 = ^(I,XiX 2 ,Xi) - (xxi,Xi +XX 2 ,X 2 + 
XY2,Y2 + XXiZ2,XiZ2) and O3 = ^(/,XiX2,/) - (xi,i + 
XXiX2,XiX2 +XXi,Xi +XX2,X2)- 

The real or imaginary counterparts of the expressions 
Eq. (fT9l) are obtained by preprocessing the noisy states 
with the corresponding toggling operations. For code C^, 
note that / and Y2 are ambiguous, and thus cannot have 
different toggler signs. On the other hand, we want them 
both to have different toggler signs than X1X2 and A'iZ2, 
which are also ambiguous. Thus, one required toggling 
operation would be: 

T+ = ^ (Hci + XiRc^Xi + YiUc^Yi + ZiUc^Zi) 

+ (XiX 2 Uc^XiX 2 + X2nciX2 

v 2 

+Y1X2RC1Y1X2 + Z1X2RC1Z1X2), ( 42 ) 

and similarly for the codes (j G { 2 , 3 }). The expres¬ 
sions obtained by pre-processing the noisy ASCs with 
unitary and toggling are 

^'(/, A 1 A 2 ,/) = Xl,I + XXiX2,XlX2 +XY2,Y2 + XXiZ2,XiZ2 
+Re{I,XiX2) + lm{Y2,XiZ2), 
ai,XiX2,Xi) =xxi,Xi +Xx2,X2 +Xy2,Y2 + XxiZ2,XiZ2 
+Re{Xi,X2) +lm{Y2,XiZ2), 
^'(I,XiX2,I) = XIJ + XXiX2,XiX2 +XXi,Xi +XX2,X2 
+ReiI,XiX2)+ReiXi,X2). 

( 43 ) 

Solving the above set of equations we have the real or 
imaginary parts of the off-diagonal terms of the process 
matrix that wre undetermined by Eq. m with out tog¬ 
gling: 

Re{I,XiX2) = ^(0'i-0'2 + 0',) = ^, 

Re(Ai, A 2 ) = ^(-O'l + O' + O') = ^, 
Im(y 2 ,XiZ 2 ) = i( 0 ( + 0 ' -0') = ^, ( 44 ) 

where O} = ^'{I,XiX2,I)-{xij+XXiX2,XiX2+XY2,Y2 + 
XXiZ2,XiZ2) ^'2 = C{ItXiX2,Xi) - {xxi,Xi +XX2,X2 + 
XY2,Y2 + XXiZ2,XlZ2) C>3 = ^'(/, A 1 A 2 ,/) - (xij + 

XXiX2,XiX2 +XXi,Xi +XX2,X2)- 
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VII. DISCUSSION AND CONCLUSION 


We developed the concept of ambiguous stabilizer 
codes, which exploit the stabilizer formalism for quantum 
error characterization rather than for quantum error cor¬ 
rection. We presented different procedures for construct¬ 
ing an [[n, fc]] ASC that ambiguously detects arbitrary er¬ 
rors on m known qubit coordinates (m < n). The Pauli 
operator basis for this set of errors forms a group. The 
ASC can be characterized as a quotient group Pm/®, 
where ® is the set of m-qubit Pauli errors ambiguous 
with the no-error syndrome. The cosets of ® form other 
ambiguous sets of errors. 

ASCs cannot be used for quantum error correction, ex¬ 
cept if the basis elements of the noise is known to have at 
most a single element in each of the ambiguous sets of the 
ASC. Quite generally, a suitable collection of ASCs can 
be employed for characterizing noise, and this is the chief 


application of ASCs. The code length for an ASC can 
be smaller than demanded by the requirement of error 
correction, making state preparations potentially simpler 
from an experimental perspective than for the techniques 
of Refs. @,0|- We developed a protocol, “quantum ASC- 
based characterization of dynamics” (QASCD), for this 
purpose, which, in comparison with the use of conven¬ 
tional stabilizer codes for CQD Q, requires smaller code 
length, but at the cost of more number of quantum oper¬ 
ations and classical post-processing. We illustrated our 
method using an example of characterization of a toy 2- 
qubit noise using three 4-qubit ASCs. 

ACKNOWLEDGMENTS 

OS and RS acknowledge financial support through 
DST, Govt. of India, for the project SR/S2/LOP- 
02/2012. OS acknowledges academic support from the 
Manipal University graduate program. 


[1] E. Knill and R. Laflamme, “A theory of quantum error- 
correcting codes,” Phys. Rev. A 55, 900-911 (1997), 
arXiv:quant-ph/9604034. 

[2] D. Gottesman, ArXiv:0904.2557. 

[3] S. Omkar, R. Srikanth, and S. Banerjee, “Characteriza¬ 
tion of quantum dynamics using quantum error correc¬ 
tion,” Phys. Rev. A 91, 012324 (2015). 

[4] M. Nielsen and I. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press 
(Cambridge), 2000). 

[5] Giacomo Mauro D’Ariano, “Quantum tomog¬ 
raphy: General theory and new experiments,” 

F. der Physik 48, 579-588 (2000) 

[6] J. B. Altepeter, D. Branning, E. Jeffrey, T. C. Wei, P. G. 
Kwiat, R. T. Thew, J. L. O’Brien, M. A.Nielsen, and 
A. G. White, “Ancilla-assisted quantum process tomog¬ 
raphy,” Phys. Rev. Lett. 90, 193601 (2003) 

[7] M. Mohseni and D. A. Lidar, “Direct 
characterization of quantum dynamics,” 
Phys. Rev. Lett. 97, 170501 (2006); “Direct char¬ 


acterization of quantum dynamics: General theory,” 
Phys. Rev. A 75, 062331 (2007) 

[8] Joseph Emerson et al, “Symmetrized char¬ 

acterization of noisy quantum processes,” 
Science 317, 1893-1896 (2007) 

[9] M. Silva, E. Magesan, D. W. Kribs, and J. Emerson, 
“Scalable protocol for identification of correctable codes,” 
Phys. Rev. A 78, 012347 (2008) 

[10] Ariel Bendersky, Fernando Pastawski, and 

Juan Pablo Paz, “Selective and efficient estimation 
of parameters for quantum process tomography,” 
Phys. Rev. Lett. 100, 190403 (2008) 

[11] Austin G. Fowler, D. Sank, J. Kelly, R. Barends, 
and John M. Martinis, ArXiv: 1405.1454; J. Gombes, 
G. Ferrie, C. Gesare, M. Tiersch, G. J. Milburn, H. J. 
Briegel, and G. M. Caves, ArXiv: 1405.5656; Y. Fuji- 
wara, ArXiv: 1405.6267. 

[12] Raymond Laflamme, Cesar Miquel, Juan Pablo Paz, and 
Wojciech Hubert Zurek, “Perfect quantum error correct¬ 
ing code,” Phys. Rev. Lett. 77, 198-201 (1996) 




